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We solve Einstein's constraint equations in the conformal thin-sandwich decomposition to model 
thin shells of non-interacting particles in circular orbit about a non-rotating black hole. We use these 
simple models to explore the effects of some of the freely specifiable quantities in this decomposition 
on the physical content of the solutions. Specifically, we adopt either maximal slicing or Kerr-Schild 
slicing, and make different choices for the value of the lapse on the black hole horizon. For one 
particular choice of these quantities the resulting equations can be solved analytically; for all others 
we construct numerical solutions. We find that these different choices have no effect on our solutions 
when they are expressed in terms of gauge-invariant quantities. 

PACS numbers: 



I. INTRODUCTION 

A 3-1-1 decomposition of Einstein's equations results in 
a set of constraint equations, which constrain the gravita- 
tional fields at all instants of coordinate time, and a set of 
evolution equations, which propagate the fields forward 
in time (e.g. [ll, l2, @]). The four constraint equations can 
constrain only a subset of the gravitational fields. There- 
fore, the constraint equations can be solved, for example 
for the construction of initial data, only after the con- 
strained variables have been separated from freely speci- 
fiable ones, and after suitable choices have been made for 
the latter (see e.g. [3, Q for reviews). 

The constrained variables are separated from the freely 
specifiable ones by choosing a decomposition of the con- 
straint equations. The conformal thin-sandwich decom- 
position iQ, |7| has been particularly popular for the con- 
struction of quasiequilibrium data; it has been used ex- 
tensively, for example, to model compact binaries con- 
taining black holes or neutron stars (see, e.g., [3, Q for 
reviews). In the conformal thin-sandwich formalism, the 
spatial metric is conformally decomposed into a confor- 
mal factor and the conformally related metric, and the 
extrinsic curvature into its trace and a traceless part. In 
the so-called extended version [7], the freely specifiable 
variables are the conformally related metric and the trace 
of the extrinsic curvature together with their time deriva- 
tives (which we may set to zero to construct equilibrium 
data), and the constrained variables are the lapse, the 
shift, and the conformal factor. 

Black holes may be constructed within the conformal 
thin-sandwich formalism by excising the black hole inte- 
rior, and imposing suitable inner boundary conditions. 
In particular, these boundary conditions may be cho- 
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sen so that the black hole is momentarily isolated, or 
in equilibrium (see [lO, llll, [l3] , also compare the isolated 
horizon formalism laid out in |13l . Il4l. 1X511 and references 
therein). As discussed in detail in [llj, these geomet- 
ric conditions lead to boundary conditions on some of 
the constrained variables in the thin-sandwich formal- 
ism, namely the conformal factor and the shift vector. 
The boundary condition for the lapse, however, remains 
arbitrary. 

Some of the choices in this formalism will clearly have 
an effect on the physical content of the solution. We can 
expect to find equilibrium solutions only if we set the 
time derivatives of the conformally related metric and 
the trace of the extrinsic curvature to zero. Also, a con- 
formally flat solution is physically distinct from solutions 
that are not conformally flat. The choice of the trace 
of the extrinsic curvature, or mean curvature, is usually 
associated with an initial temporal gauge, and the lapse 
plays a similar role. It is less clear, then, whether or how 
the mean curvature and the boundary condition on the 
lapse affect the solutions. 

In flJl, the authors found that sequences of binary 
black holes, and in particular their innermost stable cir- 
cular orbit, do depend on the horizon lapse for their ex- 
ample of a non-maximal slicing, i.e. non-zero mean cur- 
vature. This finding, however, may be an artifact of their 
particular choice of the mean curvature, namely a super- 
position of two copies of its analytical value for a single 
Schwarzschild black hole expressed in Kerr-Schild coor- 
dinates (see (I13bp below) , one for each companion in the 
binary. As the authors caution, the resulting background 
geometry then depends on binary separation, making the 
physical meaning of these sequences somewhat arguable. 

We consider a very simple physical system in order to 
analyze whether, at least in this context, the choice of the 
mean curvature and the horizon lapse affect the physical 
content of the solutions. Speciflcally, we solve the con- 
straint equations in the thin-sandwich decomposition to 
construct thin shells of non-interacting, isotropic parti- 



cles in circular orbit about a Schwarzschild black hole 
(compare [la] , whose results we generalize to account for 
the black hole). For one particular choice of the mean 
curvature and the horizon lapse the equations can be 
solved analytically (see Appendix |^ , and we construct 
numerical solutions for many others. We find that these 
different choices have no effect on our solutions when they 
are expressed in terms of gauge-invariant quantities. 



II. BASIC EQUATIONS 

A. Constraint equations 

We write the spacetime metric gat in the form 

gabdx^dx^ = -a^df + -i^j{dx' + P'dt){dx^ + f3^dt), (1) 

where a is the lapse function, /3* the shift vector, and 7ij 
the spatial metric. We further decompose the latter as 



7jj =-0 7y, 



(2) 



where ip is a conformal factor and jij a conformally re- 
lated metric. We then solve Einstein's constraint equa- 
tions in the conformal thin-sandwich decomposition (see 
[a, 0|, as well as [1, S 9 for reviews). Specifically, the 
Hamiltonian constraint becomes 



D^i; = \i}R + -^i^^K^ - -ip~'^A,jA'^ - 27rV'Vw- (3) 
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Here pN — Tab'n'^n^ is the energy density as measured by 
a normal observer, D^ — ^^^ DiDj, and Di and R are the 
covariant derivative and the Ricci scalar associated with 
the metric 7ij . We have also split the extrinsic curvature 
Kij into its trace K and a traceless part Aij according 
to 

Kij = Ay + -j,jK = ip~^A,j + -^li]K. (4) 

From the evolution equation for the spatial metric we can 
express A^^ as 



A'j 



25 ((i« 



— U'' 



(5) 



Here a = ip ^a and u^^ — dtjij, and the conformal 
Killing operator L is defined as 

{Lpy = &I3^ + &(3' - ^fWkf3^ (6) 

The momentum constraint can now be written as 
(Al^)' = {LfSy^ Dj In {a) + aD J {a'^u'^) 

+ \ai>'^D'K + 167raV/°f , (7) 

where (A^/?)* = Dj{Li3y^ is a vector Laplacian, and j* = 
—j^°-n''Tab is the mass current as measured by a normal 



observer. Finally, the trace of the evolution equation for 
Kji, combined with the Hamiltonian constraint, results 



D- 



(aV) = a^(^V~'A,,A'^ + A^4^2_^l^ (g) 
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+2^V (P + 25) - Tp^dtK + rp'D^K, 



where S = j^^Tij the trace of the spatial stress. 

The above equations form a set of equations for the 
lapse a, the shift /9' and the conformal factor tp. Be- 
fore these equations can be solved, however, we have 
to make choices for the freely specifiable quantities 7^, 
Uij = dtjij, K and dfK. For the construction of 
quasiequilibrium data it is natural to choose Uy = and 
dtK = 0. We will also restrict our analysis to spheri- 
cal symmetry, where we may assume conformal flatness, 
7jj = Vij J without loss of generality. Here 77^ is the flat 
metric in whatever coordinate system. We will, however, 
experiment with different choices for K (see equations 
P^ below), as well as with different boundary condi- 
tions for the lapse a (see (|17cp below). 

With these choices, and in spherical symmetry, the 
above equations simplify dramatically. We write the spa- 
tial metric as 

j.jdx'dx^ = V^ {dr^ + r^{dO^ + sin^ dd(f>^)) , (9) 

where r is the isotropic radial coordinate. The shift (3^ is 
now purely radial, and we abbreviate (3 = (3^ . We may 
evaluate © to find 

/ -2r2 \ 



sin^ei 



so that Aj-jA^^ becomes 



i.A- = A,2 (g^^ 



(11) 



The Hamiltonian constraint ([3]) can then be written as 



rd^ip + 2dr'ip + r 



ip^ (drP-P/ry 
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K' 



= -27rV'>Ar, 
the momentum constraint ([7]) as 

\r a ip J \ r 

= adrK+l2TnP^aj\ 

and the lapse equation ([5]) as 

7V'^ /o n /3n2 5 



(12a) 



(12b) 



,,4^2 



-h27rV'^ (p + 25) + i^^l3drK\ . (12c) 



Here we have expressed all quantities in terms of those 
variables that are used in our code. 

In the above equations the trace of the extrinsic cur- 
vature K can still be chosen arbitrarily. Following [Ifj 
we consider two different possibilities, namely maximal 
slicing 



and Kerr-Schild slicing 

2A/BH A , 2A/BH 



K 



KS 



i?2 



1 



R 



-3/2 



3Mbh \ 
R ) 



(13a) 



(13b) 



Here R is the areal radius, R = ^jPr, and we identify 
Mbh with the black hole's irreducible mass (see ([50)1 be- 
low). Kerr-Schild coordinates are identical to ingoing 
Eddington-Finkelstein coordinates. 

Throughout this paper, we use the subscript bh to refer 
to the black hole, and sh to indicate a property of the 
shell. 



B. Boundary conditions 

At spatial infinity we impose asymptotic fiatness, 
which results in the boundary conditions 



^ 



1, 



1, 



/3^0 



(14) 



as r ^ oo. 

We excise the black hole interior inside an isotropic ra- 
dius tbh and impose the black hole equilibrium boundary 
conditions of Cook and Pfeiffer [ll| on the resulting ex- 
cision surface S (compare the notion of isolated horizons 
laid out in [13, [S])- In particular, the condition 



{DiSj - Ki- 



0, 



(15) 



where s* is the outward-pointing unit normal to the hori- 
zon, ensures that this surface corresponds to a marginally 
trapped surface (apparent horizon), while the condition 



f3^\^ = (3\3r\, 



(16) 



ensures that the coordinate system tracks the horizon. 
(The tangential components of the shift vanish identically 
in spherical symmetry.) In our case, P^ becomes 



2r 6 



and (US 



M_PIl^K]] =0 (17a) 



P^\^a\s. (17b) 

As discussed in [ll|, the boundary condition for the 
lapse a is arbitrary. We will experiment with the Dirich- 
let boundary condition 



a\. 



oah, 



(17c) 



and will compare results for values of oah ranging from 
zero to unity in increments of 0.1. 



C. Matter equations 

We consider a spherically symmetric shell of isotropic, 
non-interacting particles in circular orbit about the black 
hole (compare [16|). The rest energy (baryon mass) of a 
spherically symmetric matter source may be written as 



Msu = / pou'y/^d^x = 47r / poWi^'^r^dr, (18) 



where W = 



is the Lorentz factor between 



a normal observer n"^ and an observer comoving with the 
fluid M°. For an infinitesimally thin shell we can then 
identify the rest energy density (baryon density) as 



Po = 



rSir-rsn), 



(19) 



where rgn is the (isotropic) radius of the shell. 

Since the particles are non-interacting, their stress- 
energy tensor is that of a pressureless fluid (dust): Tat = 
poUaUb- The matter sources pN, S and j^ in equations 
p2ap . (J12bp and (|12cp can therefore be expressed as 

PN = PqW^, 
S = po{W^~l), 



(20) 



The delta function in these matter sources leads to a dis- 
continuity in the derivatives of the solutions. We can 
find the jump in these derivatives by integrating equa- 
tions (O, ^ and dH) from tsh — e to rgn + e, which, in 
the limit e — > 0, results in the jump conditions 

SiyMsH^ 



drl3+ — drP- — 
dr^p+ — drip- = 

dr (oiTp): — dr (aVO- 



-MmW 
2riHV ' 



(21) 



aip 



MsH {3W^ - 2) 



Since the particles are non-interacting, their 4- velocity 
M° obeys the geodesic equation 

du^/dr + Vi^u'' = 0, (22) 

Assuming circular orbits with u^ — and du^ /dr = 
and, without loss of generality, focussing on a particle in 
the equatorial plane, we find 



r:,{uf+T^,,{u^y 



0. 



(23) 



Using the normalization condition Uqu" = — 1, which is 
equivalent to 1 = W"^ — jijU^u^ , we evaluate (|23|) to find 
the geodesic condition 



l~W'^ 



13^ Tp^P {ipdrP + 2l3drTp) - adrU \ ip'^r 



{2rdrtlj + ip) ip^ 



lA^-2. 



(24) 



for circular orbits. To compute the Christoffel symbols 
we averaged the derivatives of the gravitational field vari- 
ables inside and outside the shell (compare jlQ]). 



D. Diagnostics 

We compute the ADM and Komar masses using the 
expressions (see e.g. Chap. 7 of [17]) 



27r7aE 



(25) 



and 



Mk^^ f {D,a ~ (3^K,j) dS' , (26) 

47r J as 

where dS^ is the outward pointing unit surface element of 
a closed surface at infinity. In our case, these expressions 
reduce to 



Madm = lim — 2r dr4' 



and 



M^ 



lim r drd . 

r — >oo 



(27) 



(28) 



For all configurations considered in this paper, the two 
mass expressions are found to be equivalent to within the 
accuracy of our code when the geodesic equation ([M]) is 
used to force the particles into circular orbit (compare 
[16|). This is in agreement with a general theorem about 
the equality of ADM and Komar masses established by 
Beig [l^l and Ashtekar & Magnon-Ashtekar [l9| . We also 
define the binding energy as 



En = 



M 



ADM/K 



Ml 



BH 



Af; 



-1, 



SH 



(29) 



either terms of the ADM or Komar mass. Finally, we 
compute the black hole's irreducible mass from the area 
A of the apparent horizon. 



Mw 



{- 



1/2 



R 

2" 



^bhV'Ih 



(30) 



equations in two separate domains inside and outside the 
shell, each one represented by iV Chebyshev polynomials. 
The N coefficients can then be determined by evaluating 
the equation at N collocation points in each domain, and 
the jump conditions (j21[) can then be imposed exactly as 
matching conditions between the two sets of Chebyshev 
polynomials. 

Each set of Chebyshev polynomials is Tn{s) with < 
n < A^ — 1 and s E [—1,1]. We map the inner region into 
the interval [—1,1] with the transformation 



si 



2r'BH?'SH/'' - J'BH - ''SH 



''BH - rsK 



and the outer region with 



So = -2rsH/7' + l. 



(31) 



(32) 



Our computational domain therefore extends to r = oo, 
and we can evaluate the masses (|27p and (|28t exactly. 

In addition to the choices for q;ah and K, our solution 
depends on the parameters A/gn and rgn in p9p . and 
the excision radius tbh ■ To construct a solution of given 
black hole mass A^bh, we need to iterate over tbh un- 
til the resulting irreducible black hole mass ((5D|) agrees 
with the desired black hole mass A^bh to within a cer- 
tain pre-determined tolerance. For a given shell radius 
rgH a further iteration is needed to fix the Lorentz factor 
W in such a way that the solution satisfies the geodesic 
condition (|24p . and the particles are in circular orbit. In 
practice, we instead fix W and then iterate over rgn until 
([M)) is satisfied. 

For each case we start the iteration from the analytical 
solution for K = and ckah = provided in Appendix 
\X[ and continue until the solution has converged. Specif- 
ically, our convergence criterion requires that the relative 
change between iteration steps in any of the fields is less 
then 10~^° at all collocation points. 



Tests 



III. NUMERICS 



A. Code 



We developed a pseudo-spectral code to solve the dif- 
ferential equations (fT^ subject to the boundary condi- 
tions (HD) and (fT7|) as well as the jump conditions (PTI) . 
using Chebyshev polynomials as basis functions (sec [20] 
for a recent review of spectral methods) . Equations (J12bp 
and (|12cp can be solved directly, while equation (|12ap has 
to be linearized and then solved iteratively. 

One complication arises as a consequence of the jump 
conditions ([2T|) . Representing the solution functions 
across these jumps as a linear combination of the contin- 
uous Chebyshev polynomials would result in undesirable 
Gibbs phenomena. To avoid this problem, we solve the 



We tested our program for a number of different known 
vacuum solutions expressing the Schwarzschild geometry 
in different coordinate systems, as well as an analytical 
solution describing thin shells around black holes that we 
derive in Appendix |^ 

As a first vacuum test (for which we set Afgn = and 
rgH = lOAfBH), we considered Schwarzschild in "stan- 
dard" isotropic coordinates, representing the symmetry 
plane in a Carter-Penrose diagram. In our code, we can 
produce the solution 

M M 

^ = 1 + —, a^^l-—, 13^0 (33) 
2r 2r 

by choosing K — Kus — and a ah = 0. After applying 
the transformations (PT|) and ([5^ . these solutions be- 
come linear in our code's coordinates si and sq, meaning 
that the solution can be represented exactly in terms of 
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FIG. 1: The average error over all collocation points E 
as a function of the number of collocation points N for 
two analytic representations of the Schwarzschild geometry 
in isotropic coordinates. The label AP denotes the analytic 
puncture solution presented in [2111 . while the label KS de- 
notes the Kerr- Schild solution in [ll[. The errors drop off ex- 
ponentially as a function of the number of collocation points 
A'^, until a 'floor' of specified tolerance has been reached. 



FIG. 2: The conformal factor ip and the lapse a as a function 
of areal radius for the analytic shell solution {K = and 
QfEH ~ 0) for a Lorentz factor oi W — 1.20 and a mass ratio 
of Mbh/Msh = 1, corresponding to an areal shell radius of 
about 8.013AfBH. The field variables are continuous across 
the shell, but, according to the jump conditions (|21|l their 
derivatives are not. Our numerical solutions agree with the 
analytical ones to within better than 10"^", making the error 
far smaller than the line width in the graph. 



the first two Chebyshev polynomials. For any A^ > 1, our 
code therefore converges to the correct solution vifithin 
the predetermined tolerances. 

More interesting arc the two isotropic representations 
of the Schwarzschild geometry presented in [11[ and 
[2l| . The former is a transformation of Kerr-Schild 
(Eddington-Finkelstein) coordinates to isotropic coordi- 
nates, keeping the same time slicing, which we can pro- 
duce by choosing K = Kks and a ah = l/'x/Z in our 
code. The latter is an isotropic representation of a max- 
imal slice (with the critical parameter C — 3a/3M^/4, 
see [23, u3|), which has recently attracted interest as 

I r I ij / 5 J 

an analytic "puncture" solution (compare [2j|). We can 
produce this solution by choosing K — Kus = and 
ttAH = 3v3/16 in our code. 

To measure the deviation from an analytic solution, 
we compute the average of the absolute error e„ at all N 
collocation points in each of the two domains, 



j^N ^ 



1 



2N 



2N ^-^ 

n=l 



(34) 



where the collocation points TV -I- 1 < n < 2N are in 
the region outside of the shell. In Fig. [T] we graph this 
error as a function of the number of collocation points 
A^. As expected, the errors fall off exponentially for all 
variables, until they reach a floor corresponding to the 
predetermined tolerance. 

In addition to these vacuum solutions we also con- 



sider an analytic solution describing thin shells of non- 
interacting particles around a static black hole. As we 
demonstrate in Appendix [^ we can solve the differential 
equations P^ subject to the boundary conditions (fTH) 
and (|17p as well as the jump conditions (pij) analytically 
for maximal slicing and a ah = 0. As for the solution 
(j33p . the field variables become linear in our code's vari- 
ables Si and sq (see (|A2p and (jA3p below), so that they 
can be represented exactly for any A'^ > 1. In addition 
to testing the solution of the field equations, however, 
this test also verifies that our code correctly solves the 
jump conditions at the shell. As in example, we show in 
Fig. [5] the analytic and numerical solutions for the lapse 
a and the conformal factor -ip as a function of areal ra- 
dius R for a Lorentz factor ofW — 1.20 and a mass ratio 
Mbh/Msh = 1- Our numerical solutions agree with the 
analytical ones to within better than 10""'^'', making them 
indistinguishable in the plot. 



IV. RESULTS 

We now construct constant-mass sequences, mean- 
ing sequences of varying shell radius rgn but constant 
shell rest mass Msh and black hole irreducible mass 
A/bh- Specifically, we focus on "extreme-mass-ratio" 
sequences with Mbh/Msh — 1000 and equal-mass se- 
quences Men/AfsH = 1- For both choices of the mass 
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FIG. 3: The ADM binding energy ([29} Eb as a function of 
the shell's area! radius R for an extreme-mass-ratio sequence 
with Mbh/Msh ~ 1000. Values of the horizon lapse oah 
range from zero to unity in increments of 0.1. The graph 
includes eleven solid lines (representing an evaluation of the 
binding energy (|29p in maximal slicing) and eleven dashed 
lines (for Kerr-Schild slicing) . All 22 lines coincide within our 
numerical error. We computed these sequences using N = 26 
coUocation points in each domain. 



ratio we construct sequences for our two choices of the 
extrinsic curvature (|13[) . maximal slicing and Kcrr-Schild 
slicing, and for the horizon lapse (|17cp ranging from 
ttAH = to QfAH = 1 in increments of 0.1. 

In the following, we will graph the binding energy (j29p 
as a function of the areal radius. Typically, the bind- 
ing energy Mb is only a small fraction of the involved 
masses, meaning that the relative error in the binding 
energy is larger than that for the masses - and hence the 
fields - themselves. We found that, to achieve similar ac- 
curacy, the extreme-mass-ratio sequences (for which the 
binding energy is much smaller than for the equal-mass 
sequences) required slightly more collocation points than 
the equal-mass sequences. 

In Fig.[3]we show the binding energy for extreme-mass- 
ratio sequences with A/bh/^^sh = 1000. The graph rep- 
resents 22 plots, corresponding to the eleven different val- 
ues of aAH and to evaluating the ADM binding energy 
(P5)l for both maximal slicing and Kerr-Schild slicing. To 
within the accuracy of our code, all 22 lines agree with 
each other, so that they all lie on top of each other and 
appear as one line in Fig. [31 

The minimum of the binding energy corresponds to the 
innermost stable circular orbit (ISCO). In the extreme- 
mass-ratio limit we may neglect the particles' self-gravity, 
so that we are effectively solving for a test-particle in 
the Schwarzschild geometry. As expected, we find that 
the ISCO is located at i? = 6Mbh, representing another 



FIG. 4: The ADM binding energy Eb as a function of 
the shell's areal radius R for an equal mass sequence, 
Mbh/Msh ~ 1- As in Fig. [S] values of a ah again range 
from zero to unity in increments of 0.1, and the graph again 
contains eleven lines each for Maximal and Kerr-Schild slic- 
ing. All 22 lines again coincide within our numerical error. 
We obtained these results with A'^ = 20 collocation points in 
each domain. 



independent test of our code. 

In Fig. [H we show the equivalent graphs for equal-mass 
sequences with Mbh/A^sh = 1. As for the extreme-mass- 
ratio sequences in Fig. [3] all 22 graphs coincide to within 
our numerical accuracy. We note that the ISCO is now 
located at a larger radius of about 9.367Mbh- 

Even in this case, we do not find any evidence that the 
choice of the slicing condition p3)) . or the choice of the 
boundary condition for the lapse on the horizon (|17cp . 
have any effect on the physical content of our solutions. 
Clearly, these different choices lead to different solutions 
for the conformal factor ip, the lapse a and the shift (3. 
When expressed in terms of gauge-invariant quantities, 
however, all our solutions become indistinguishable to 
within the accuracy of our numerical code. 



V. SUMMARY 

We solve Einstein's constraint equations in the ex- 
tended conformal thin-sandwich decomposition to con- 
struct spherical shells of non-interacting, isotropic par- 
ticles in circular orbit about a non-rotating black hole. 
We construct these solutions for both maximal slicing 
and Kerr-Schild slicing (see equations p3|) ). and for a 
number of different choices for the horizon lapse. These 
different choices lead to very different solutions for the 
lapse a, the shift /?' and the conformal factor ip. How- 
ever, when expressed in terms of gauge-invariant quanti- 



ties - for example the binding energy as a function of the 
shell's areal radius for given shell and black hole masses 
- our solutions become indistinguishable. At least in the 
limited context of our spherically symmetric solutions, 
these findings provide no evidence that the choices for 
the mean curvature and horizon lapse affect the physical 
content of the solutions. 
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APPENDIX A: AN ANALYTICAL SOLUTION 

For maximal slicing {K — 0) and the lapse bound- 
ary condition aAH = 0, we can find an analytic solution 
to the differential equations (fT^ . subject to the bound- 
ary conditions (fT4)l and (flT)) as well as the jump condi- 
tions pT|) . The solution depends on the input parame- 
ters MsH, W, tbh, and rsn- We may then enforce the 
geodesic condition (|M|) to eliminate one of these vari- 
ables, and compute the mass of the black hole from ([30)1 . 
The solution presented here represents a generalization 
of the solutions of 16], who considered the same system 
but without a black hole. 

We begin with the momentum constraint (|12bp . For 
K = Q and q:ah = we find that 



/3 = 



(Al) 



is a self-consistent solution to both the equation and its 
boundary conditions. This implies that Eqs. (|12ap and 
(I12c[) for respectively the conformal factor ij) and the com- 
bination ail) reduce to flat Laplace equations in the vac- 
uum regions away from the shell. In spherical symmetry, 
the only possible solutions are of the form ki -\- ^2/?", 
where fci and k2 are arbitrary constants that have to be 
determined from the boundary conditions. For each func- 
tion we need four conditions to determine these constants 
both in the interior and the exterior of the shell. These 
four conditions arise from the outer boundary conditions 
(|14p . the inner boundary conditions (J17p . continuity of 
the functions at the shell, and the jump conditions (PT|) 
for their first derivatives. Using these conditions, we find 



(' + ^)' 



r < rsH 



iP= { , , A2 

1 H , r>rsii 



for the conformal factor and 
c{l-rBH/r) 



a ■ 



a{l + rsii/r) ' 

1 + (c (rBH - rsii) - rsii) /r 



r < rsH 



r > rsH 



(A3) 



1 + (a {rsu + '"bh) - rsn) /r 
for the lapse. Here the constants a and c are given by 

MshW 



a== 1 



2»'shV'sh 



and 



c= 1- 



MsH^SH {SW^ - 2) 
rsH (iWips-a) 



(A4) 



(A5) 



Inserting these constants, which themselves depend on 
the values of the conformal factor and the lapse at the 
shell, into (U2]) and (TOl). we find 



; 1 M , ^BH 

2 V ^SH 



and 



2MsiiW 

TBH + 7-SH 



1/2- 



asH = P 



{3+p)W^-2 
2rs-iiips-iiW 



(A6) 



(A7) 



Here we have abbreviated p = (rsn + 7'bh)/(?'sh — tbh)- 
So far these solutions depend on all four parameters 
-MsH, W, tbh, and rsH. We now find a relation between 
these parameters by inserting (jA2p and (jA3P into the 
geodesic equation (^51) . which yields 



r!^- 



(4)pr- 



(A8) 



_ ^ShC (rBH - t-sh) (t'sh (g - c) + tbh {a + c + 2ac)) 
2(a(rsH + ?'bh)) (a^BH - ^sh) 

Unfortunately, this expression is not very useful for our 
purposes in this form. We find an alternative form by 
evaluating ([M)) for our solution (|A2p and (jA3p . which 
results in a fifth order polynomial for W, 



= 



4(rBH - ?'sh)^?'sh 

+Msh(?'bh - rsB)^{rBii + rsB)W 

-2rsH(?'BH - '■SH)(5r|H - 14rBH?'SH + 5rgjj)W^^ 

-2MsH(rBH - 2rsH)(?'BH - ?'sh)(?'bh + rsB)W^ 

+6rsH(rBH - ^sh)(?'bh - 4rBH?'SH + rlaW^ 

+MsH(rBH - 2rsH)'(»^BH + rsii)W^ (A9) 



In the limit of a zero-mass black hole we have tbh = 
and (IA9I1 reduces to 



-4rsH + MsbW + lOrsHM^^ - AMshW^ 
-6rsH W^"^ + 4MsHTy ^ , (AlO) 
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which agrees with the corresponding equation (26) of [Ta| ■ 
Instead of parameterizing the solution by tbhi it is 
more desirable to fix the black hole mass Mbh = Mi„- 
Towards this end we combine ((X9|) with ((X2|) and ([5U)) . 
which results in a quadratic equation for rsH- The solu- 
tions to this equation are 



^-SH 



4MBHr|H + 2MBHMsHrBHW - 2M^^rBH 
2(M|h - 2MBHr-BH)' 



(All) 



We henceforth ignore the "— " solution, as only the "+" 
solution refers to stable solutions. 

Finally, we may insert (|A11[) back into (|A9p . which 
results into a polynomial for tbh that can be factored 
into two cubic polynomials 







(asrliii 



3 
BH 



a2r|H 



ai^^BH +ao) X 



(bsr 



^ar-BH + hrBH 



(A12) 



Here we have abbreviated tbh = ^ehZ-^sh, and the 
coefficients ai and bi are, in terms of the mass ratio 
q = Mbh/Msh, 



0.3 

ai 



b3^32q^{-2 + 3WY 
-mq{-2 + 3VF2)[-V + GqW 
+ (6g2 - 1)W^ - 9qW^ + 2W^] 
-4g(-2 + 3VF2)2(8q2 - 8qW + 5W^) 
2[16q'^ - ASq^W + (4V " 'iSq^)W'^ 
+ {lUq^ - l2q)W^ + (36/ - 136g2 
(38g - 108q^)W^ 



(A13) 



6 _ . 

6i = 2{-2 + 3W^)[-8q'^ 



[lOdq^ 



2)W^ - SOqW'^ + W^] 



-{12q^ - 20q^)W^ 



- IGq^W 
{lOq ~ 24:q^)W^ 
+ {26q^ - 2)W'^ - 12qW^ + 3W^] 
ao = bo^-qW^{-2q + W + 3qW'^ -2W^f 



We can now construct a solution for given masses Mbh 
and MgH and Lorentz factor W as follows. Given the 
mass ratio q and W we first find the six solutions for 
tbh from the two polynomials in (jA12[) . We then insert 
the corresponding values tbh hito (jAlip . which yields 
six solutions for rgfj- The largest real solution is the 
solution of interest; we keep only this solution as well as 
its corresponding value of tbh and disregard all others. 
These values can then be inserted into (|A2[) and (IA3p . 
which determines the solution completely. 
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